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SYMBOLS 


[A] = 

[M] = 

d_ 

dt 

t 

t » «u» rv> p 

AE 

R 

T 

T 0 »T = 

X A 

F' 

F 


OJ 




T 


tj(T) = 
k. 

i 

f. 


Concentration of specie "A" 

"Third body" concentration 

Rate of change with respect to time 

Time 

= Exponents signifying order of reaction 

Activation energy of reaction 

Universal gas constant 

Temperature 

Mean temperature 

Mole fraction of specie "A" 

Pre- exponential factor in Arrhenius reaction rate expression 

Same as F T , but associated with rate expression given in terms of 
mole fractions 

Frequency of sinusoidal oscillation 
Amplitude of sinusoidal oscillation 
Period of sinusoidal oscillation 
Rate of reaction at temperature T 

Rate of reaction associated with free radical reaction for specie "i" 
Rate of photolysis associated with specie "j" 
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EFFECT OF TEMPERATURE OSCILLATION 
ON CHEMICAL REACTION RATES IN THE ATMOSPHERE 


INTRODUCTION 

There has in the recent decade been a very considerable interest in the chemis- 
try of the earth's atmosphere, and many photochemical models have been pub- 
lished. Rather than review the field of photochemical modelling, I refer the 
interested reader to a recent review article by H. U. Diitsch (1973), which gives 
a good introduction to the photochemistry of atmospheric ozone and also has a 
bibliography containing the names of most of the scientists active in the field. 

In a very interesting study (1973) Johnston and Whitten examined the production 
and destruction of atmospheric ozone in terms of the Chapman photochemical 
theory, ozone transport, and reactions with oxides of nitrogen. They calculate 
that the globally averaged ozone production rate by the Chapman mechanism is 
five times larger than the globally averaged destruction rate by the Chapman 
mechanism. Johnston and Whitten (1973) also calculate ozone transport to the 
troposphere as being about 1% of the gross rate of ozone formation, using num- 
bers for transport given by Brewer and Wilson (1968), 

Johnston and Whitten (1973) show that the excess ozone production rate by the 
Chapman mechanism can be readily balanced if destruction of ozone by oxides 
of nitrogen is introduced. A modest background concentration of nitric oxide is 
assumed to exist in the atmosphere. The explanation is certainly plausible but 
its validity must await confirmation until sufficiently accurate and extensive 
measurements of atmospheric nitric oxide become available. 

The purpose of the present paper is to consider the possible effect of tempera- 
ture oscillation on effective reaction rates. The nonlinearity of the Arrhenius 
expression for the chemical rate is such that the mean rate under a sinusoidal 
temperature fluctuation will be larger than the rate at the mean temperature. 

Let us examine how large the difference is likely to be. The author came across 
the problem some ten years ago when the possible effects of turbulent tempera- 
ture and density fluctuations on chemical reaction rate measurements needed to 
be estimated (Glassman and Eberstein, 1963). Now, the effect of temperature 
fluctuations on atmospheric ozone chemistry will be examined. 


ANALYSIS 

The chemical reaction for specie "A" is given by an expression of the form 
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( 1 ) 


— [A] = 

d t 


F'[A] g [B]™ tC] n exp 



where 


p - £ + m + n 


( 2 ) 


is the order of reaction. Any of the exponents may be zero. The factor F may 
be a weak function of temperature, but is generally taken to be constant. 


In the atmosphere, one generally deals with constant pressure levels for which 
total density is inversely proportional to temperature. Equation (1) may thus 

be written: 


d [A] • 
dt 







(3) 


Let us now subject the system described by Eq. (3) to a sinusoidm tempera -ire 
oscillation, and determine the ratio of "mean rate" to "rate at mean tempera- 
ture.” The mole fractions X, and the quantity "F" will be assumed to be unaf- 
fected by the temperature fluctuations. We thus have: 


T 

— = 1 + /3 sin <*>t , 


(4) 


where fi is oscillation amplitude, and: 


V (T) 

V CO 


({’ 


- A E 

T _f> [1 + /3 sin a>t] p exp PrY 


f-AE 1 T r 

|_RT 0 1 + fi sin a>tj 


dt 


(5) 


| Tq P exp -r^] 


dt 


The quantity r/ (T) is the temperature sensitive part of 


d [A] 
dt 


i.e. 
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„ F-ae" 

7,(T) = T p exp |~rj~ 


( 6 ) 


V (T) is the rate at the mean temperature T, and 77 (T) isjdie mean rate integrated 
over a cycle for the fluctuating temperature T. T 0 and T represent the same 
quantity. 'Equation (5) was evaluated numerically for a cycle and some of the 
results are shown in Figures (1) and ( 2 ). In Figure 1 the ratio rfi T)/?? (T) is 
plotted versus the fluctuation amplitude for several values of the overall reaction 
order P and the quantity A E / R T. 

Figure 1 shows considerable enhancement of the mean reaction rate for tem- 
perature fluctuations exceeding 10%. 

Figure 2 shows the reaction rate ratio versus AE/ET 0 for P = 0 for tempera- 
ture fluctuations of 10%, 15%, 20% and 25%. It is seen that significant reaction 
rate enhancement occurs as AE/RT 0 goes above 10. 

Let us now examine what temperature oscillations and activation energies one 
may expect in the atmosphere. Most upper atmospheric reactions involve either 
photochemical dissociation or free radical recombination. These processes 
have a very weak temperature dependence. 

However, reactions like 


0 3 + 0 - 20 2 


do have anon-trivial temperature dependence, i.e. (Dutsch, 1973) 


k = 


6 x 10' 1 2 exp 



giving a value of approximately 10 to the quantity AE/R T 0 . 

In examining the Chapman mechanism, Johnston and Whitten (1973) give the fol- 
lowing expression for the steady state concentration which ozone is potentially 
capable of attaining: 
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Figure 1. Ratio of Mean Rate to Rate at Mean Temperature 
Versus Temperature Fluctuation Amplitude 
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V(T> 
V (T) 



Figure 2. Ratio of Mean Rate to Rate at Mean Temperature Versus Activation Energy 



( 7 ) 


Ms = /* 2*2 CM1V 

[0 2 ] V f 3 k 3 / 

with 

.4 / S02\ 

k 2 = 1. 10 x 10' 34 exp 


and 


1.33 x 10- 11 



The following reactions have been considered: 

f 2 

0 2 4- hv - 0 + 0 

k 2 

M + O + 0 2 -» 0 3 + M 

f 3 

0 3 + hv - 0 2 + 0 

k 3 

0 + °3 °2 + °2 


(1) 

( 2 ) 

( 3 ) 

(4) 


At 50 km it is reasonable to take the mean temperature as 260°K, and the fluc- 
tuation as ±10% (Theon et al., 1972). 

The ratio of photochemical rates associated with the fluxes f 2 and f 3 is assumed 
to remain unaffected by the oscillation in Eq. (7). The formation rate associated 
with k 2 [M] is increased by 3% and the destruction rate associated with k 3 is 
increased by 8.5%. Referring to the quantity in brackets in Eq. (7), one would 
expect it to decrease by 5.5% resulting in a 3% drop of steady state ozone. 

Referring to the table given by Johnston and Whitten (1973), the following values 
are obtained from the Chapman mechanism: 
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Global rate of ozone production = 500 
Global transport to troposphere = 6 
North-South transport ± 4 

The above values are in units of 10 29 molecules/ second. 

Now* a 3% change in the global production rate is 15 units, or more than the total 
North-South transport, and more than twice the global transport to the tropo- 
sphere. 

We may conclude that the effect of temperature oscillation on atmospheric 
photochemistry is comparable to transport effects, but much less than the factor 
of 5 required by Johnston and Whitten (1973) to balance the ozone produced in 
the Chapman mechanism. 

The above analysis is in many ways preliminary in nature, and a more thorough 
study is underway. However, the author felt that it is of great importance to 
make atmospheric modellers aware of the fact that chemical reaction rates in 
the atmosphere respond in a nonlinear and non-negligible fashion to temperature 
oscillations in the atmosphere. 
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